Laser-induced nonsequential double ionization at and above the 
recollision-excitation-tunneling threshold 
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We perform a detailed analysis of the recoUision-excitation-tunneling (REST) mechanism in laser- 
induced nonsequential double ionization (NSDI), in which the first electron, upon return, promotes 
a second electron to an excited state, from which it subsequently tunnels, based on the strong-field 
approximation. We show that the shapes of the electron momentum distributions carry information 
about the bound-state with which the first electron collides, the bound state to which the second 
electron is excited, and the type of electron-electron interaction. Furthermore, one may define a 
driving-field intensity threshold for the RESI physical mechanism. At the threshold, the kinetic 
energy of the first electron, upon return, is just sufficient to excite the second electron. We compute 
the distributions for helium and argon in the threshold and above-threshold intensity regime. In the 
latter case, we relate our findings to existing experiments. The electron-momentum distributions 
encountered are symmetric with respect to all quadrants of the plane spanned by the momentum 
components parallel to the laser-field polarization, instead of concentrating on only the second and 
fourth quadrants. 



I. INTRODUCTION 

Electron-electron correlation in strong laser fields has 
raised considerable interest for over a decade, in particu- 
lar in the context of laser-induced nonsequential double 
and multiple ionization Jj. Concrete examples are the 
early measurements of a "knee" in the double ionization 
yield as a function of the laser-field intensity, which devi- 
ates from the predictions of sequential models in orders of 
magnitude 01, and the peaks in the electron momentum 
distributions in nonsequential double ionization (NSDI), 
as functions of the electron components Pn\\ [n = 1,2) 
parallel to the laser-field polarization Such peaks 
occur at nonvanishing parallel momenta and cannot be 
explained by a sequential mechanism. 

Presently, it is an established fact that NSDI occurs 
due to the inelastic recoUision of an electron with its par- 
ent ion [4] . In this recollision, the first electron gives part 
of the kinetic energy it acquired from the driving field to 
a second electron, which is then freed. 

The simplest type of recollision which can lead to this 
phenomenon is electron-impact ionization. Thereby, the 
first electron, upon return, provides the second electron 
with enough energy so that it is able to overcome the 
second ionization potential of the target in question and 
reach the continuum. Both electrons leave simultane- 
ously and lead to distributions peaked at nonvanishing 
momenta, occupying the first and third quadrant of the 
plane Pi||P2|| spanned by the parallel momentum compo- 
nents. 

Electron-impact ionization is possibly the most exten- 
sively investigated NSDI mechanism [5-[l3|. This is due 
to the fact that, until the past few years, it was suf- 
ficient to describe key features observed experimentally 
in the electron momentum distributions. Concrete ex- 
amples are the peaks near the non-vanishing momenta 



Pi\\ — P2\\ — i2-^/E^, where Up is the ponderomotive en- 
ergy, and the recently reported V-shaped structure 
which is a signature of the long-ran ge c haracter of the 
electron-electron interaction [tI IsI. HoI Il2l ITsj. Moreover, 
from the theoretical viewpoint, this mechanism is consid- 
erably easier to model, as compared to other rescattering 
processes. This is particularly true in the context of semi- 
analytical models, such as the strong-field approximation 
[10,1 El, El- 
Recent experimental results, however, reveal other 
rescattering mechanisms apart from electron-impact ion- 
ization. For instance, if the driving-field intensity is be- 
low the threshold intensity, the kinetic energy of the 
recoUiding electron is no longer sufficient to make the 
second electron overcome the ionization potential of the 
singly ionized core and reach the continuum fl^ EE] • In 
this intensity range, the energy transferred to the core 
is only enough to raise the second electron to an excited 
bound state, from which it subsequently tunnels. This 
process is known as recollision-excitation-tunneling ion- 
ization (RESI) . This mechanism is also important for spe- 
cific targets, such as argon [17|], even if intensities above 
the threshold are taken, provided the driving pulse is long 
enough. 

For this specific mechanism, the first electron leaves 
near a crossing of the driving field, whilst the second 
electron is freed near a field maximum. The time delay 
between both electrons suggests that they leave with op- 
posite momenta. This implies that the second and the 
fourth quadrant of the Pi||P2|| plane are expected to be 
populated . A population in such regions of the paral- 
lel momentum plane has also been observed for NSDI of 
diatomic molecules, especially if the molecules in ques- 
tion are aligned perpendicular to the laser-field polar- 
ization [19l - l22| . Furthermore, it has been recently shown 
that electron-impact ionization, alone, would not be very 
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useful for retrieving the molecular structure in the experi- 
mentally relevant intensity range [23| . All these examples 
suggest that the RESI mechanism is becoming increas- 
ingly important, in view of the driving field intensities 
and targets employed. 

This mechanism, however, is considerably less under- 
stood than electron-impact ionization. In fact, apart 
from early results in which only the RESI yield has been 
calculated within a semi-analytical framework [isj , most 
existing results are the outcome of computations, either 
classical [16, 19, 20, 24, 25] or quantum mechanical [23 |. 
for which the different rescattering mechanisms cannot 
be easily disentangled. 

In previous work , we have studied the RESI phys- 
ical mechanism, within the strong-field approximation. 
In this framework the transition amplitude is written as 
a multiple integral, with a time-dependent action and 
slowly-varying prefactors. This approach has the partic- 
ular advantage of providing a clear space-time picture of 
the process in question, in terms of electron trajectories, 
while retaining features such as quantum-interference ef- 
fects (27,] . 

We have shown that the RESI mechanism may be un- 
derstood as the combination of two processes. For the 
first electron, a behavior similar to rescattered above- 
threshold ionization is present. The main difference lies 
on the fact that part of the kinetic energy of the first elec- 
tron is given to the parent ion. This leads to a maximal 
kinetic energy slightly smaller than the high-order ATI 
cutoff lOUp. The second electron, which is tunnel ion- 
ized at a subsequent time, behaves in the same way as 
in direct above-threshold ionization. Hence, its maximal 
kinetic energy is given by the direct ATI cutoff, i.e., 2Up. 
This allowed one to determine kinematic constraints for 
the RESI mechanism, which lead to distributions which 
equally occupy the four quadrants of the Pi\\P2\\ plane. 
For an early derivation of somewhat different kinematic 
constraints for this mechanism see, e.g., Ref. [2^]. 

In the above-stated investigations, however, we as- 
sumed the prefactors in the strong-field approximation 
amplitude to be constant. Nonetheless, in a more re- 
alistic scenario, one should take into account that such 
prefactors cause a bias in momentum space. This may 
lead to electron momentum distributions concentrated on 
the second and fourth quadrants of the Pi||P2|| plane. 

This work is organized as follows: In Sec. |lll we will 
briefly recall the SEA transition amplitude for the RESI 
mechanism, which has been derived in [2^ . We will start 
from the general expressions (Sec. Ill Ap . which will be 
solved by saddle-point methods. The pertaining saddle- 
point equations are derived in Sec. Ill Bl Subsequently, 
in Sec. Ill CI we will provide the specific prefactors em- 
ployed for the hydrogenic systems to be investigated in 
this work. In Sec. IIIII we employ this approach to com- 
pute electron momentum distributions for Helium and 
Argon. For the latter species, an explicit comparison 
with the results in Ref. [14| is performed. Finally, in Sec. 
IIVI we state the main conclusions of this paper. 



II. TRANSITION AMPLITUDE 

A. General expressions 

The transition amplitude describing the recollision- 
excitation-tunneling ionization (RESI) mechanism, 
within the SEA, can be written as (for details on the 
derivation see |26j). 
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Eq. ([T]) describes the physical process in which an 

electron, initially in a bound state \ijjg'^^ >, is released 
by tunneling ionization at a time t" into a Volkov state 
\k{t') >. Subsequently, this electron propagates in the 
continuum from t" to a later time t' . At this time, it 
is driven back by the field and rescatters with its parent 
ion. In this collision, it excites the second electron, which 

(2) (2) 

is bound at \tpg >, to the state \^pe >, through the 

interaction V12. Finally, the second electron, which is in 

(2) 

a bound excited state |V'e >, is released by tunneling 
ionization at a time t into a Volkov state |p2 (t) >. The 
final electron momenta are described by p„(n — 1,2). 
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In the above-stated equations, Eng {n = 1,2) give the 
ionization potentials of the ground state, E^f. (n = 1,2) 
give the ionization potentials of the excited state and 
V{ri) and Vion(''2) correspond to the atomic binding 
potential of the system as seen by the first and second 
electron, respectively. 

The form factors ([3]) and (O contain all the informa- 
tion about the binding potential of the first and second 
electron, respectively. The form factor (|4]) contains all 
the information about the interaction of the first elec- 
tron with the singly ionized atom, which is an inelastic 
interaction. If the electron interaction is only dependent 
on the difference between both electron coordinates, i.e., 
if Vi2(ri,r2) = Vi2{vi — r2), Eq. (HJ may be rewritten as 

_ li2(Pi-k) 

Vp.eMa - (2^)3/2 

X j d3^2e-^(P^-'')-'-n^f^(r2)]*Vf ^(r2),(6) 

with 

Fi2(pi-k) = ^^-^ j d3rexpH(pi-k).r]yi2(r) (7) 

and r = ri — r2. 

Clearly, Vkg and Ip^e are gauge dependent. In fact, in 
the length gauge p„ (r) = Pn -I- A(t) and k(T) = k + 
A{t){t = t', t"), while in the velocity gauge p„ (r) = p„ 
and k(r) = k. This is a direct consequence of the fact 
that the gauge transformation xi^v = exp[— iA(T) • r] 
from the length to the velocity gauge causes a translation 
p — >■ p — A(t) in momentum space. Due to the fact that 
these shifts cancel out in Eq. ([S]), Vpie,kg remains the 
same in the length and velocity gauges. 

B. Saddle-point analysis 

The multiple integral in Eq. ([1]) will be solved us- 
ing saddle-point methods (for details see Ref. p9|). 
For that purpose, we must find the coordinates 
{ts,t'^,t",'ks) for which ^(pi, P2, k, t, t', t") is stationary, 
i.e., for which the conditions dtS{pi,p2,^,t,t' ,t") = 
dt'Sipi,p2,k,t,t',t") = dt"S{pi,p2,k,t,t',t") = and 
i9k'S'(Pi, P2, k, t, t' , t") ~ are satisfied. This leads to the 
equations 

[k-f A(0]'--2i?i„ (8) 

[pi + A{t')f = [k + A{t')f - 2{E2g - E2e). (10) 
and 

[P2 + A{t)]^ = ~2E2e (11) 



Eq. ([8]) gives the conservation of energy at the time 
i", which, physically, corresponds to tunneling of the first 
electron. Since tunneling has no classical counterpart, 
this equation possesses no real solution. Eq. © con- 
strains the intermediate momentum k of the first elec- 
tron, so that it can return to the parent ion. Eq. (jlOp ex- 
presses the fact that the first electron returns to its parent 
ion at a time t' and rescatters inelastically with it, giving 
part of its kinetic energy Eret{t') = [k + A(t')]^ /2 to the 
core. Under this interaction the second electron is excited 
from a state with energy E2g to a state with energy E2e ■ 
The first electron leaves immediately and reaches the de- 
tector with momentum pi. Finally, Eq. ([TT|) describes 
the fact that the second electron tunnels at a later time 
t from an excited state of energy E2e- 

1. Momentum constraints 

The saddle-point equations (ITUl) and pD) provide use- 
ful information on the momentum-space regions popu- 
lated by the RESI mechanism, and on the shapes of the 
electron-momentum distributions. For instance, Eq. (fTTj) 
is identical to that describing tunnel ionization for direct 
above-threshold ionization. This implies that the maxi- 
mal kinetic energy of the second electron at the detector, 
if the field can be approximated by a monochromatic 
wave, is roughly given by 2Up. 

Furthermore, the electron is leaving the excited state 
with largest probability when the electric field E{t) = 
—dA{t)/dt is maximum. If the time dependence of the 
laser field is such that A{t) vanishes when E{t) is at its 
peak (for instance, monochromatic fields), then 

-2y%<p2<2y%. (12) 

If, to first approximation, we neglect the momentum 
components perpendicular to the laser-field polarization, 
one can see that the momentum of the second electron, in 
the parallel momentum plane, is expected to be centered 
around vanishing momentum p2\\ and be limited by the 
bounds p2\\ — ±2-y/C^. 

One should note that, due to the fact that Eq. pT|) de- 
scribes a tunneling process, there is no classically allowed 
region for the second electron. A non-vanishing perpen- 
dicular momentum, effectively, will lead to an increase in 
the potential barrier and a suppression in the yield. In 
fact, Eq. (fTTj) can also be written as 

[P2\l + A{t)]^ ^ ~2E2e, (13) 

where E2e = E2e—P2± is an effective ionization potential. 

The saddle-point equation (fTO| . on the other hand, 
yields information about the momentum of the second 
electron. According to this equation, 

- A{t) - ^2Edis <Pi\\< -A{t) + ^2Eiis, (14) 

where E^m = E^,,{t' ,t") - {E2g - E2e) - plJ2 and 
-£'kin(^', t") denotes the kinetic energy of the first electron 
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upon return. For a monochromatic field, the electron 
returns most probably near a crossing of the laser field, 
one may use the approximation A{t) ~ I^JlTp in the 
above-stated equation. In this case, we also know that 



(max) 



the kinetic energy £'kin(i',i") < 3.17C/p. Hence, iS^jfj" 



< 



/2 and 



< -2a 



Eq. (1151) allows one to delimit a region in momentum 
space for piy centered around — 2-y/T7p and bounded by 

'2E^^^\ In contrast to the previous case, there may 
be a classically allowed region for the momentum of the 
first electron if the parameters inside the square root 
are positive, i.e., if 3.17Up > (i?2g — E2e) + Pij_/2. 
For increasing perpendicular momentum and/or bound- 
state energy difference, this region will become more and 
more localized around —2y/lJp until it collapses. There- 
fore, it is also possible to distinguish a threshold and 
an above-threshold behavior in the context of recoUision- 
excitation-tunneling. One should note, however, that in- 
tensities below the recoUision-excitation threshold (i?2g — 
E2e) — 3.l7Up do not make physically sense, as the en- 
ergy of the returning electron would not be sufficient 
to promote the bound electron to an excited state. If 
(i?2g - E2e) < 3.17C/p the well-known cutoff of lOUp for 
rescattered above-threshold ionization is recovered [3ij |. 

In view of the above-mentioned constraints, the ex- 
pected maxima of the electron momentum distribution 
are located at the most probable momenta {pi\\,P2\\) = 
{±2,yUp, 0), and, after symmetrizing with respect to the 

exchange pi ■(-> P2, at {pi\\,P2\\) = (0, ±2y/C7p). This im- 
plies that, if the field can be approximately described by 
a monochromatic wave, the outcome of our model should 
be distributions in the Pi\\P2\\ plane, which are symmet- 
ric upon p„ — — p„, n = 1,2 and upon pi o P2, and 
which equally occupy the four quadrants of the parallel 
momentum plane. 



2. Bound-state singularity 

Finally, due to the saddle-point equations ^ and 
(ITT]) , for exponentially decaying bound states the pref- 
actors (|3]) and ([5]) exhibit singularities in the length- 
gauge formulation of the strong-field approximation. 
This is due to the fact that these prefactors will be 

inversely proportional to (jk + A{t")]^ + 2Eig^ and 

([p2 -I- A(t)]2 -t- 2i?2e)'" , where n, m are integers. For the 
problem addressed in this specific work, however, only 
the prefactor Vp^e will influence the shape of the elec- 
tron momentum distributions. The prefactor Vkg will 
affect the electron momentum distributions only quanti- 
tatively. Hence, to first approximation, one can consider 
Eq. ([3]) as constant. A similar problem for the electron- 
impact ionization mechanism in NDSI has been discussed 
in detail in 0]. 



To overcome the singularity in Vp^e, one needs to em- 
bed this prefactor into the action, which now reads 

^(pi , P2, k, i, t', t") = 5(pi, P2 , k, t, t', t") - I In Vp,e. 

(16) 

This will lead to modifications in the saddle-point equa- 
tion 9t5(pi, P2, k, t, t' ,t"), which is now given by 



[P2 + A{t)]^ = -2E2e + idt In Vp, 



(17) 



The main consequence of such a modification is that the 
drift velocity of the second electron is no longer pure 
imaginary. This will lead to a splitting in the ionization 
time t for each orbit, as compared to the non-modified 
case. Depending on the velocity in question, the barrier 
the electron must tunnel through in order to reach the 
continuum will either widen or narrow. This means that, 
with regard to the non-modified action, ln\[t] will either 
increase or decrease. 



C. Prefactors 

In this work, we are particularly interested in expo- 
nentially decaying, hydrogenic bound states. This means 
that, in general, the bound-state wavefunction reads 



(18) 



where n, I and m denote the principal, orbital and mag- 
netic quantum numbers, the index a refers to the electron 
in question, and the angular coordinates are given by 9a 
and if a- In this case, the binding potentials V(ri) and 
Vion{r2) wiU be given by 



(19) 



where Va yields either V or yion, and Z^q corresponds 
to the effective electronic charge. The general expres- 
sions for the prefactors in this work are provided in the 
appendix. 

Below, we state the specific prefactors to be employed 
in Sec. IIIIl for Helium and Argon. In the former case, 
upon collision, the second electron may be excited from 
the Is state to either the 2s or the 2p state, while in the 
latter species it may undergo a transition from the 3p 
state to the 4s or the 4p state. One should note that the 
prefactor V^^e is gauge dependent. In the length gauge, 
P2(t) = P2 + A(t), and P2{t)\\ = P2\\ + M^) while, in 
the velocity gauge, P2(i) = P2 and P2[t)\\ — P2\\- The 
prefactor Vpie.kgj on the other hand, is gauge invariant. 



1. Excitation Is — >■ 2s 

Let us first consider the simplest case, in which the 
second electron is excited to 2s. This gives the prefactors 
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with 



{K^,E2g,E2e) = l^{^j2ETg+2^2E2e) + {2E., 



2g) 



n3/2 



-2(2£;2e)'/'-6£;2ey2^. (22) 



and 



aE2g,E2e)^V^e+V2E, 



2S- 



(23) 



The above-stated equations can also be written in 
terms of the momentum components parallel and per- 
pendicular to the laser field polarization, denoted hy Pa\\ 
and Pa±, {a — 1,2), respectively. In this case, 



T/(2s) 
P2e 



[[P2(t)|| 



pl^ + 2E2eV 



(24) 



{k-Pl\\Y + pi^,E2g,E2e 



[{k-p,\^Y+pl^+e{E2g,E2eW 

(25) 



2. Excitation Is — >■ 2p 

If, on the other hand, the second electron is excited 
to 2p, one must consider three degenerate states, corre- 
sponding to the magnetic quantum numbers m — ±1,0. 
This yields 



T/(2p) 
P2e 



Mt)Y 



and 



y{lS^2p) 

Pie, kg 

with 



, 2-[n"(%.,¥'pJ]* (26) 

i2E2e + [p2{tf) 

Vi2{pi - k)7?2(A^^ E2g, E2e) [Fr(0K, f.)]*, 

C{E2g,E2e)V^ 



?72(k ,E2g,E2e) 



iCHE2g,E2e) + n^y 



(27) 
(28) 



Since the electron may be excited to any of the 2p states, 
we will consider the coherent superposition 
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V2py 



^2p, 



(29) 



where (r2 
plies that 



(2) \ _ ,/,(2) 
2p, 



V'y, (i"2), with j = x,y,z. This im- 



T/(2p) 
P2e 



[P2(t)] 



2£;2 



[P2(i)]^ 



r/3(P2(t)) 



(30) 



(ls->2p) 



Pie, kg 



Fl2(pi - k)772 (K', £;2g, £^2e)/3(«:), (31) 



where the angular dependency is given by 

/3(q) = (sin(?q cos(^q + sin6'g sincpq + cos^g). (32) 

Thereby, we employed the usual relations between spher- 
ical polar coordinates and the spherical harmonics. 

One may write the above-stated expressions in terms 
of the electron momentum components parallel and per- 
pendicular to the laser-field polarization. In this case, 
Eq. (Uni reads 



T/(2p) 
P2e 



V2{t)\ 



/3(p2(t)). (33) 



(2E2e+[p2{t)l\] 

In /?(p2(i)), the angles Op^ and ipp^ are given by 

P2it)\\/^J[p2it)\\f +PI 



arccos 



(34) 



and (/Pp2 — SLTCCOs[p2{t)x/p2it)±], respectively. In Eq. 
(pit , n"^ = {k — Pi\\Y + Pi_L and the angles 6*^ and (^^ 
read 



arccos 



k~Pl\\) /-y/(fc -pi||)2 -f p2_L 



(35) 



and = arccos[pi2;/pij_], respectively. This angular 
dependence will be washed out when the transverse mo- 
mentum components are integrated over (see Sec. IIIip . 



3. Excitation 3p — > 4s and 3p — 4p 

Finally, we will assume that the second electron, ini- 
tially in 3p, will be excited either to the 4s or to the Ap 
state. Similarly to the procedure adopted in the previ- 
ous section, we will consider a coherent superposition of 
the 3pa;, Zpy and ipz states for the initial state of the 
electron, i.e., 



V3 ( 



(36) 



If the electron is excited to the 4s state, the excitation 
prefactor V^^^~^g'^^ will exhibit an angular dependence 

given by /3(k), and the tunneling prefactor Vp'^l' will not 
depend on the angular variables. Both prefactors also 
have a radial dependence on P2{t) or k. If, however, the 
electron is excited to the 4p state, one must take the final 
state as 



4? 



1 



73 ( 



^4p„ 



*S0) 



(37) 
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i.e., as a coherent superposition of Ap^, Apy and Apz. In 
this case, the angular dependence of Vp^e will be embed- 
ded in /3(p2(t)). The angular dependence oiV^^~^g^^ will 
be more complex and will involve the sum of the orbital 
angular momenta of the two electronic bound states in- 
volved. Due to the higher quantum numbers involved, 
the prefactors are messier than those in the previous sec- 
tions and will not be written down explicitly. They can, 
however, be obtained from the general expressions in the 
appendix. 

III. ELECTRON MOMENTUM 
DISTRIBUTIONS 

In this section, we will compute electron momentum 
distributions, as functions of the momentum components 
(Pi II 7^211) parallel to the laser-field polarization. We ap- 
proximate the external laser field by a monochromatic 
wave, i.e., 

E(i) So sin wtGa;. (38) 

This is a reasonable approximation for pulses whose du- 
ration is of the order of ten cycles or longer (see, e.g. Q 
for a more detailed discussion). In this case, the elec- 
tron momentum distributions, when integrated over the 
transverse momentum components, read 

J^(Pi||,P2||) = yy"dV±rfW|Mfl(Pi,P2) (39) 

-f AfL(pi,P2) +Pl ^ P2P, 

where Affl(pi,P2) is given by Eq. ([!]) and d?Pni_ = 
Pn±dpn±d(pp^ . If A{t ± T/2) = -A{t), where T = 27r/a; 
denotes a field cycle, the actions Sl, Sr corresponding 
to the transition amplitudes Ml and A/^ obey the sym- 
metry 5L(pi,p2,i,t',<") - 5i?X-Pi,-P2,t±r/2,t'± 
T/2,t" ± T/2). The distributions have also been sym- 
metrized with respect to the exchange pi -O- P2. To 
a good approximation, the quantum-interference terms 
[My(pi,P2)]*M^(pi,P2),J^ 7^ /i, get washed out upon 
the transverse-momentum integration, so that it is suf- 
ficient to add the above-stated amplitudes incoherently. 
Here, we considered V^g as constant and we integrate the 
transition amplitude over the azimuthal angles Lpp^ . 

We will now briefly discuss how the prefactors Vp^e 
and Vpje,kg behave with regard to the integration over 
(fp^ . Obviously, if the second electron is excited from an 
s state to an s state, the prefactors Vp^e and Vp^e.kg do 
not depend on this parameter. However, if a transition 
from or to a p state is considered there will be an angular 
dependence in such prefactors. 

For instance, tunneling ionization from a p state would 
lead to the argument /3(p2(t)) in Eq. ([50]) . If excitation 
from an s state to a p state or vice versa takes place, the 
angular dependence of the prefactor V^je.kg is given by 
P{k) in Eq. (|3T|) ). When integrated over the azimuthal 
angles (pp^ , |/3(k)| and |/3(p2)| will yield 27r, so that the 



angular dependence of these prefactors can be neglected. 
If, however, states of higher orbital quantum numbers 
are involved, or if the initial and excited bound states of 
the second electron are p states, this dependence will be 
more complex. 

The simplest scenario is if all prefactors are nonsin- 
gular, such as in the velocity- gauge formulation of the 
SEA. In this case, they will contribute to the electron 
momentum distributions as |Vp2e|^ and iVp^e.kgP- In the 
length-gauge SEA, however, Vp^e exhibits a singularity, 
and therefore must be incorporated in the action accord- 
ing to Eq. ([T7)) . This singularity, however, is present 
only in the radial part of this prefactor. Therefore, the 
angular parts are still slowly varying, and may be treated 
as above. The radial part of the prefactor, however, must 
be incorporated in the action. 



A. Intensity dependence 

We will commence by having a closer look at how 
the momentum-space constraints affect the electron mo- 
mentum distributions for different driving-field intensi- 
ties. Eor that purpose, we will assume that the prefac- 
tors Ipickg and Vp^e are constant, and vary the laser- 
field intensity. Eor the lowest intensity, the kinetic en- 
ergy of the returning electron is just enough to pro- 
mote the second electron to an excited state, i.e., we 
are considering the recollision excitation (RESI) thresh- 
old E2g — — 3.17J7p. This intensity, however, is be- 
low the electron- impact ionization threshold, i.e., > 
3.l7Up. The intermediate intensity has been chosen such 
that i?2g — < 3.17?7p, i.e., above the threshold for 
recoUision-excitation. Nevertheless, this intensity is not 
sufficient to make the second electron overcome the ion- 
ization potential and be freed by electron-impact ioniza- 
tion. Einally, the highest driving-field intensity consid- 
ered in this section is far above the recollision-excitation 
threshold, and slightly above the electron-impact ioniza- 
tion threshold. This implies that rescattering is classi- 
cally allowed for both physical mechanisms. The compu- 
tations in this section have been performed for Helium, 
and the pertaining results are presented in Eig. [TJ 

As an overall feature, the distributions exhibit four 
peaks at {pj\\,p^ii) = (±2-y/[7p, 0), with j,!/ — 1,2 and 
j V. These peaks agree well with the constraints dis- 
cussed in the previous sections. This holds even if the 
driving-field intensity is just enough to excite the sec- 
ond electron and the only allowed momenta are ±2y^I^ 
[Fig. [11(a)]. Physically, this means that the first elec- 
tron will reach its parent ion most probably at a crossing 
of the driving field, reaching the detector with the most 
probable momenta of ±2^/17^, while the second electron 
will reach it with vanishing momentum. 

The shapes of the distributions, however, differ con- 
siderably. Indeed, at the RESI threshold intensity 
[Fig. [11(a)], one observes ring-shaped distributions. As 
the intensity increases, this distributions become more 
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FIG. 1: Electron momentum distributions for Helium (Eig = 
0.97 a.u., E2g = 2 a.u. and i?2e = 0.5 a.u.) in a linearly 
polarized, monochromatic field of frequency uj = 0.057 a.u.. 
In the picture, we considered all prefactors to be constant. 
Panels (a), (b) and (c) correspond to a driving- field intensity 
I = 2.16 X W/cm^ 7 = 2.5 X W/cm^ and / = 3 x W/cm^ 
respectively. The contour plots have been normalized to the 
maximum probability in each panel. 



and more elongated along the Pn\\ axis [Fig. [11(b)], un- 
til the maxima merge and cross-shaped distributions are 
observed [Fig. [TJ(c)]. 

This change of shape may be understood by analyz- 
ing the momentum-space constraints. The widths of the 
distributions are determined by the tunnel ionization of 
the second electron from an excited state. This process 
has no classical counterpart and leads to distributions 
peaked at p2\\ — and which vanish at p2|| = i2-^/C7p, 
i.e., at the direct ATI cutoff. Increasing the intensity will 
only make the effective potential barrier smaller or wider, 
and thus affect the overall yield, but will not change such 
constraints. 

The elongations in the distributions are determined by 
the rescattering of the first electron. This rescattering, 
in contrast, delimits a momentum region which is highly 
dependent on the driving-field intensity. Therefore, its 
width in momentum space will vary. Specifically, at the 
RESI threshold, there will be maxima in the distributions 
at pi|| = ±2,yUp due to the fact that the first electron 
rescatters most probably at a field crossing. However, as 
these are the only classically allowed momenta, the dis- 
tributions will be fairly narrow around this value. With 
increasing driving-field intensity, the classically allowed 
region defined by Eq. ()15[) will become more and more 
extensive and this will cause the elongation. 

Note that the electrons are indistinguishable so that 
the above-stated arguments hold upon the exchange 
Pi\\ P2\\- Hence, the horizontal and vertical axis in 
the parallel momentum plane will be equally affected. 



B. Bound-state signatures 

We will now investigate how the shape of the bound 
state to which the second electron is excited is imprinted 
on the electron momentum distributions. We will also 
employ different gauges and types of electron-electron 
interaction. Explicitly, we will assume that the second 
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FIG. 2: Velocity-gauge electron momentum distributions for 
Helium {Eig — 0.97 a.u., E2g = 2 a.u. and E2e = 0.5 a.u.) 
in a linearly polarized, monochromatic field of frequency uj = 
0.057 a.u. and intensity I = 2.16 x W/cm^. In panels (a) 
and (c), the first electron has been excited to 2s, while in 
panels (b), and (d) it has been excited to 2p. The interaction 
employed is indicated in the figure. The contour plots have 
been normalized to the maximum probability in each panel. 



electron is either excited by a contact-type interaction 
^12 ^(''i ^ ^2) = '^(''1 " ''s) or by a long-range. Coulomb 
type interaction V"/^^(ri — r2) = l/(ri — r2). In the for- 

while in the latter case 



mer case, V^/2^(pi — k) = const. 

y/2''^(pi-k)^l/(pi-k)2. 

In order to perform a direct comparison, we will take 
the same parameters as in Fig. [Tl but incorporate the 



prefactors V^ll'^g'^ and V^te , or V^'l'^g^' and V^^e 
corresponding to the Is — >■ 2s or Is — >■ 2p excitation with 
subsequent tunneling, respectively. 

In Fig. [21 we consider the lowest intensity in the pre- 
vious figure and the velocity gauge. If the electron 
is excited to the 2s state [Fig. [21(a)], we observe four 
spots which are slightly elongated along the p„|| axis. 
Hence, in comparison to its constant prefactor counter- 
part, i.e.,Fig. [Tl(a), there was a narrowing. This narrow- 
ing is caused by the interplay of two features in the pref- 

(2s) 

actor Vp2e • First, this prefactor exhibits two symmetric 
nodes, which, for vanishing transverse momentum are lo- 
cated at p2\\ = i^yVp. As the transverse momentum 
increases, these minima move towards vanishing paral- 

(2s) 

lei momenta. Second, Vpje decreases very steeply with 
transverse momenta P2j.. Hence, upon integration over 
this parameter, the main contributions will be caused by 
small values of p2± and will vanish near p2\\ = -ii^/Up. 
If, on the other hand, one assumes that the second 
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electron is excited to 2p, there is both a broadening in 
the distributions and a sphtting in their peaks. These 
features are depicted in Fig. [2j(b). The sphtting occurs 
at the axis Pn\\ = 0, n = 1,2, and is caused by the fact 
that Vp^e exhibits a very pronounced node at vanish- 
ing momenta, i.e., exactly where one expects Im[i] to be 
minimum and the yield to be maximum. This has been 
verified by a direct inspection of the radial dependence 
of Eq. (|30p . and omitting the Vp^e prefactor in our com- 
putations. The latter procedure caused the additional 
minima to disappear (not shown). The broadening in 
the distributions as compared to the Is —J' 2s case is a 
consequence of the much slower decrease in Vp^e^ with 
increasing transverse momentum p2± and of the absence 
of the nodes at P2\\ = ^^JTTp. There are also additional 
nodes at the diagonal = P2|| f^nd at the anti-diagonal 
Pi|| = "P2|| of the pi||P2|| plane. 

In this intensity regime, there seems to be little dif- 
ference in the shapes of the distributions if the electron 
is excited to the 2s state, regardless of whether the first 
electron interacts with its parent ion through a contact or 
a Coulomb interaction [Figs. [21(a) and (c), respectively]. 
This is possibly caused by the fact that the prefactor 
lp^2e\ due to its fast-decaying behavior, delimits a very 
narrow region in momentum space. This adds up to the 
very restrictive momentum constraints. In contrast, the 
effect of the Coulomb tail is much more critical if the elec- 
tron is promoted to the 2p state. Indeed, for a Coulomb 
type interaction [Fig. [11(d)], the splitting of the peaks 
at the axis p„|| = remain, but the nodes at px\\ = Pi\\ 
and p\\\ = ^V2\\ disappear as compared to its contact- 
interaction counterpart [Fig. [21(b)]. This is caused by 
the fact that the former minima are a characteristic of 
the prefactor, whereas the latter are mainly deter- 

mined by momentum-space effects. The Coulomb inter- 
action introduces a further momentum bias, and washes 
out the latter nodes. 

We will now discuss what happens if the intensity of 
the driving field is such that 3.17C/p > Eig — E^e- These 
results are displayed in Fig. [31 for the highest intensity 
in Fig. [1] As expected, all distributions are much more 
elongated along the axis Pn\\ — 0, as compared to the 
low-intensity case. The imprint, however, of the different 
bound states to which the second electron is excited and 
from which it subsequently tunnels are the same as in 
the below-threshold regime. Indeed, we notice that there 
is a narrowing in the distributions for the Is — 2s case 
[Figs. [31(a) and (c)], and a splitting in the peaks at the 
four axis in the Is — 2p case [Figs. [31(b) and (d)]. This is 

not surprising, as the prefactors Vpje'' and Y^^e exhibit 
the same functional dependencies as before. 

The shapes of distributions, however, change much 
more critically in this intensity regime, with regard 
to the type of electron-electron interaction, than for 
the intensity used in Fig. [2l For all cases, the dis- 
tributions computed using the Coulomb-type interac- 
tion (Figs. [31(c) and [31(d)) are much more localized in 
the low-momentum regions than those computed with a 
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FIG. 3; Velocity-gauge electron momentum distributions for 
the same parameters as in the previous figure, but driving- 
field intensity / = 3 x W/cm^. In panels (a) and (c), the first 
electron has been excited to 2s, while in panels (b), and (d) it 
has been excited to 2p. The interaction employed is indicated 
in the figure. The contour plots have been normalized to the 
maximum probability in each panel. 



contact-type interaction (see Figs. [31 (a) and [31(b)). This 
is expected, as 1^2 (pi — k) favors low momenta for the 
former, while it is constant for the latter. Physically, this 
reflects the fact that rescattering of the flrst electron is 
now allowed to occur over an extensive region in momen- 
tum space. Hence, it does make a difference whether the 
second electron is excited by a long-range or zero-range 
interaction. 

We will now perform an analysis of the electron- 
momentum distributions in the length gauge. In this 
case, the prefactor Vp^e governing the tunneling of the 
second electron exhibits a singularity, and must be incor- 
porated in the action. The modifications in the action 
read, for the 2s and 2p bound states. 



^5tlnyp(2|) 



.2£;(i)p2(i)||([P2(t)||]'+pL-6£;2 



X+(p2(i))x-(P2(0) 



and 



(40) 



P2e 



. E{t){^\^^2E2,-i{p2{M^) 

P2{t)\\ ([p2(i)||]'+pL+2i?2< 



, (41) 



respectively, with 

X±(P2(<)) = ([P2(t)||]'+pL±2£;2e) and E{t) 
-dtA{{). 
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FIG. 4: Length-gauge electron momentum distributions 
for Helium in a linearly polarized, monochromatic field. 
Throughout, we assumed Vp^^cWg ~ const., V12 to be a 
contact-type interaction, and incorporated Vp^l^ in the ac- 
tion. Panels (a) and (c) correspond to the trajectories for 
which the barrier has been narrowed by the modifications in 
the action, while Panels (b) and (d) correspond to those for 
which it has been widened. The upper and lower panels cor- 
respond to driving field intensities / = 2.16 x 10^''W/cm^ and 
/ = 3 X 10^''W/cm^, respectively. In order to perform a quan- 
titative comparison, we are providing the explicit values for 
the NSDI yield. 



For each orbit, the tunneling time of the second elec- 
tron will split into two values, as compared to the non- 
modified action. This has particularly important conse- 
quences as far as Im[t] is concerned, since it provides a 
rough measure of the width of the barrier through which 
the second electron tunnels. Physically, this means there 
will be one set of orbits for which the effective potential 
barrier will be widened, and another one for which it will 
be narrowed. 

In Fig. m we present the contributions of each of the 
orbits resulting from this splitting for a final 2s state, for 
two different driving-field intensities. For simplicity, in 
order to single out the effect of the modified action, we 
took the rescattering prefactor Ip^ckg to be constant. 

In general, the distributions differ quantitatively in a 
factor between 1.5 and 1.7, depending on whether lm[t\ 
decreased [Figs.HJa) and (c)] or increased [Figs.Sfb) and 
(d)]. This shows that the splitting in this quantity is 
small, and therefore both contributions are comparable. 

Furthermore, the distributions displayed in Fig. |4] are 
strikingly similar to those observed in Fig. [T] (see panels 
(a) and (c) therein), for which only constant prefactors 
have been considered. Indeed, the width of all distribu- 
tions, along the axis, is determined by the direct ATI 
cutoff, i.e., — 2y/c7p < pn\\ < I^JTTp. At first sight, this is 



unexpected, as we are assuming that the second electron 
is tunneling from a 2s state. As previously discussed, the 

(2s) . . . / 

prefactor Vpae exhibits a node in p2|| = ^yUp, which 
leads to a narrowing of the distributions along the p„|| 
axis. An inspection of Eq. (l24l) also suggests that, were 
it not for its singularity, the length-gauge prefactor would 
be very similar to the velocity-gauge prefactor. This is 
a consequence of the fact that the second electron is 
leaving when the field E{t) is near its maximum. For 
a monochromatic field, this implies that the vector po- 
tential A(t) is practically vanishing. 

One should note, however, that we are considering only 
the individual contributions from each of the orbits orig- 
inating from the modification of the action. It is very 
likely that, in order to recover the structure determined 
by the prefactor Vp'^e' , one must consider the coherent su- 
perposition of all the orbits originating from the splitting 
of Im[t] when computing the yield. Since these contribu- 
tions are comparable, one expects the above-mentioned 
nodes to be recovered due to quantum-interference ef- 
fects. 



C. Comparison with experiments 

We will now perform a direct comparison with the re- 
sults in Ref. [Tjl. In particular, in this reference, the 
distributions encountered have been modeled employing 
the electron-impact ionization physical mechanism and 
a modified ionization threshold for the second electron. 
Apart from that, however, in view of the driving-field in- 
tensities involved, one expects recollision-excitation tun- 
neling to be present. 

For that purpose, we will consider argon and the same 
laser-field parameters as in Ref. [l3| (c.f. Fig. 2 therein). 
We will assume, however, that, when the first electron 
recoUides, it excites the second electron from the 3p state 
either to the 4s or to the 4p state. Thereby, we took the 
velocity gauge, and assumed that the first electron inter- 
acts with the ion by a Coulomb or contact interaction. 

The results for the 3p — 4s excitation are presented 
in Fig. [51 An overall feature in the distributions are two 
main maxima along the p„||, n = 1,2, axis. These fea- 
tures are mainly caused by the Vp^e^ prefactor for the 
tunnel ionization of the second electron, which decays 
very rapidly with increasing transverse momenta and ex- 
hibit nodes near p2\\ = iO.Sy^C^. In general, we have 
verified that this prefactor determines the shape of the 
electron-momentum distributions. Secondary maxima, 
around one order of magnitude smaller, occur due to the 
rescattering prefactor Vp^l~^g^\ This prefactor exhibits 
an annular shape around pi|| = P2|| = 0- 

The existing experiments, however, do not lead to 
distributions concentrated along the axis of the Pi||P2jj 
plane. The results for Helium in the previous section 
suggest that a p state may lead to broader distributions. 
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FIG. 5: Velocity-gauge electron momentum distributions for 
Argon in a linearly polarized, monochromatic field of fre- 
quency ij — 0.057 a.u. The electron is excited from 3p to 
4s, i.e., Eig = 0.58 a.u., £29 = 1.02 a.u. and E2e = 0.40 
a.u. in our calculations. The laser-field intensity in panels 
(a) and (c), and panels (b) and (d) is I = 9 x 10"W/cm^ 
and I — 1.5 x lO^^W/cm^, respectively. The type of interac- 
tion V12 taken is indicated in the figure. The contour plots 
have been normalized to the maximum probability in each 
panel. We have verified, however, that the highest yields on 
left-hand panels are between one and a half and two orders of 
magnitude smaller than those on the right-hand side. 



FIG. 6: Velocity-gauge electron momentum distributions for 
Argon in a linearly polarized, monochromatic field of fre- 
quency iu — 0.057 a.u. The electron is excited from 3p to 
4p, i.e., Eig = 0.58 a.u., E2g = 1.02 a.u. and £2^ = 0.31 
a.u. in our calculations. The laser-field intensity in panels (a) 
and (c), and panels (b) and (d) is / = 9 x lO^^W/cm^ and 
/ — 1.5 X lO^^W/cm^, respectively. The type of interaction 
V12 taken is indicated in the figure. The contour plots have 
been normalized to the maximum probability in each panel. 
We have verified, however, that the highest yields on left-hand 
panels are between one and two orders of magnitude smaller 
than those on the right-hand side. 



For that reason, we will assume that, instead, the second 
electron is excited to the 4p state. 

Fig. |6] depicts the electron- momentum distributions for 
Argon under the assumption that the electron was ex- 
cited from 3p to 4p. All distributions in the figure exhibit 
four main maxima, which are broader than those in Fig.[S] 
and almost split at the axis Pn\\ = 0. These maxima are 

mainly determined by the prefactor Vp^eK which has a 
node at the axis for low transverse momenta and nodes 
around P211 = ± y^lTp across a wide transverse- momentum 

range. Apart from that, the prefactors V^^^~^g^^ decay 
more slowly with regard to the transverse momenta. This 
implies that, upon integration, a larger momentum region 
will be contributing to the NSDI yields. As in the previ- 
ous case, this prefactor also leads to secondary maxima 
(see Figs. EKc) and (d) for concrete examples). In all 
cases, both in Figs. [5] and El a Coulomb type interaction 
mainly introduces a bias towards lower momenta. 

Despite the above-mentioned broadening, the electron- 
momentum distributions in Fig. [B] are still considerably 
narrower than those observed in Ref. [3]. Within our 

framework, this constraint is imposed by the Vp2e^ pref- 
actor. In fact, we have verified that, for large prin- 



cipal quantum number, this prefactor always exhibits 
nodes at lower absolute momenta than the ATI cutoff 
of P2\\ = '^\fUp- In fact, if Vp^e is taken to be constant, 
the distributions become considerably broader and a bet- 
ter agreement with the experiments is obtained. This is 
shown in Fig. [71 as ring-shaped distributions with four 
symmetric maxima at pj^n = p2\\ ^^^d pj^n = — P2||- Such 

maxima are mainly determined by the V^J^e^g''' pref- 
actor. One should note, however, that this procedure 
is inconsistent from a theoretical perspective: Since the 
electron has been excited to the 4p state, it should subse- 
quently tunnel from it. Hence, the pertaining prefactor 
must be taken. 



IV. CONCLUSIONS 

Our analysis of the rescattering-excitation ionization 
(RESI) mechanism shows that the NSDI electron mo- 
mentum distributions depend on the interplay between 
the relevant momentum-space regions, the type of in- 
teraction exciting the second electron, and the spatial 
dependence of the bound states involved. We will com- 
mence by discussing each of these issues separately. 



11 



(a) 

Contact 




-4 -2 2 4 -4 -2 2 4 



FIG. 7; Velocity-gauge electron momentum distributions for 
Argon in a linearly polarized, monochromatic field of fre- 
quency Lj — 0.057 a.u. and intensity / = 1.5 x 10^''W/cm^. 
The electron is excited from 3p to 4p. We have taken the 
prefactor Vp^.e to be constant. The type of interaction V12 is 
indicated in the figure. The contour plots have been normal- 
ized to the maximum probability in each panel. 



The shapes of the electron momentum distributions, 
are determined by the interplay between two different be- 
haviors, associated with the collision of the first electron 
and the tunneling of the second electron. The momentum 
region determined by the tunnel ionization of the second 
electron from an excited state will be always restricted 
by the direct ATI cutoff. The relevant momentum region 
will not change regardless of the driving-field intensity, 
as this will always be a classically forbidden process. 

The first electron, on the other hand, rescatters inelas- 
tically with its parent ion, giving part of its kinetic energy 
upon return to excite the second electron. Hence, if its 
maximum return energy is larger than the energy differ- 
ence i?2g — E2e, rescattering has a classical counterpart. 
This implies that there will be a classically allowed re- 
gion in momentum space. If, however, this energy is just 
enough to excite the second electron, the classical region 
will collapse. Hence, the extension of the relevant region 
in momentum space related to the rescattering of the 
first electron will depend on the driving-field intensity. 
Hence, the distributions become increasingly elongated 
as the intensity increases. 

This also implies that one may define a threshold 
driving-field intensity for the RESI mechanism. This in- 
tensity is considerably lower than that necessary for the 
second ionization potential to be overcome by the second 
electron, i.e., for electron-impact ionization to occur. 

Apart from that, we have observed that the bound 
states involved in the process leave very distinct finger- 
prints on the electron momentum distributions. This is 
particularly true for the bound state of the second elec- 
tron, prior and subsequently to excitation. In fact, the 
widths of the distributions, their shapes and the number 
of maxima present will strongly depend on the princi- 
pal and orbital quantum numbers of the bound states 
involved. 

In contrast, the type of interaction V12 by which the 



second electron is excited influences such distributions in 
a less drastic way. Indeed, a long-range. Coulomb inter- 
action mainly introduces a bias towards lower momenta, 
as compared to a contact-type interaction. 

A very important observation is that all distributions 
encountered in this work are equally spread over the four 
quadrants of the Pi||P2|| plane. Under no circumstances 
have we found electron momentum distributions concen- 
trated only on the second and fourth quadrant of this 
plane, as reported in the literature [la. [20I [21|. 

Within our framework, the above-stated symmetry can 
immediately be inferred from Eq. (|40|. Nonetheless, 
one could argue that our approach does not include the 
residual binding potential in the electron propagation in 
the continuum. Recent results, however, from a classical- 
trajectory computation in which the Coulomb potential 
has been incorporated, also revealed the same symmetry 
if only the RESI mechanism is singled out [25*1. This is 
a strong hint that our results are not an artifact of the 
strong-field approximation. 

Hence, we suspect that, in the existing literature, the 
contributions from the RESI mechanism to nonsequential 
double ionization also equally occupy the four quadrants 
of the Pi||P2|| plane. They may, however, be difficult to 
extract, as explained below. 

In many situations addressed in the literature, the 
driving-field intensity is high enough for electron-impact 
ionization to occur. This means that this latter NSDI 
mechanism is also present, and fills the first and third 
quadrant of the Pi||P2|| plane. Since, in many ab ini- 
tio models, the different rescattering mechanisms are 
difficult to disentangle, the contributions from electron- 
impact ionization possibly obscure those from RESI in 
this region. In the second and fourth quadrant of the 
parallel momentum plane, the former contributions are 
absent and those from RESI can be more easily identi- 
fied. In our approach, electron-impact ionization is ab- 
sent from the start. 

If the driving-field intensities are below the electron- 
impact ionization threshold, the second electron may 
no longer be provided with enough energy to overcome 
the second ionization potential. Consequently, RESI 
becomes more prominent and the distributions equally 
occupy the four quadrants of the parallel momentum 
plane. In fact, ring-shaped distributions centered around 
Pi II — P2II =0 have been observed experimentally for this 
intensity region [l3, [l^ ■ 

Our results are far more localized near the Pn\\ — 
axis than the experimental findings. This discrepancy 
may be due to the following reasons. First, for higher 
intensities employed in Ref. [l3|, collisional excitation 
may take place not only to the 4s or to the 4p state, but 
also to highly lying states, or to a coherent superposition 
of excited states. To take this into account may be needed 
in order to reproduce the experimental data. 

Second, at the relevant driving-field intensities, one ex- 
pects the excited states to be distorted by the field, and 
the propagation of the electron in the continuum near 
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the core to be influenced by the residual ionic potential. 
This implies that a semi-analytical treatment beyond the 
strong-field approximation is necessary (see, e.g., [l3l.l30j 
for other phenomena and the electron-impact ionization 
case, respectively). Such a treatment is outside the scope 
of the present paper, and will be the topic of future work. 
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V. APPENDIX 

In this appendix, we provide the general expressions 
for all prefactors employed in this paper. We will make 
no simplifying assumption on the initial states of the first 
and second electron, and on the excited state to which 
the second electron is promoted, apart from the fact that 
they are given by hydrogenic wave functions. In order to 
compute the prefactors, we will employ the expansion 



1=0 m=-l 

Yr{o,.,^,.) (42) 

where q denotes a generic momentum, Tq the coordinate 
of the a*^ electron, and the spherical Bessel func- 
tions of the first kind. This expression will be both used 
in the derivation of Vp^e and Up^e^kg, together with the 
orthogonality relation 



i^r(0„,(p„)J Yr{ec,.vc.)dn = 5iv5„,^., (43) 

where denotes the solid angle. 

For the former prefactor, Eq. ([S]) reduces to 

Vp,e - JJ eM-iiP2{t)-^2)]Rniir2)Yr{02,ip2)r2dr2dn 
= 4^(-z)'rr(<?p.,^pJIi, (44) 

nOO 

Ii^ r2Rni{r2)ji{p2{t)r2)dr2. (45) 



with 



Similarly, Eq. ^ reads 

Vpie.kg-Ui2(pi-k)l2, (46) 

with 

I2 = yrf=^r2e*(''-Pi)-^i?„,,,(r2) [y™=(02,^2)]* 

i?n,;,(r2)r"'(^2,¥'2), (47) 



where the indices g and e in the principal, orbital and 
magnetic quantum numbers refer to the ground and ex- 
cited states, respectively. 

We will now compute the radial integrals Ii and I2 
explicitly. For that purpose, let us consider a generic 
Hydrogenic radial wavefunction 



Rniir) = Cnir'eM-y/2E^r] 



^ 'lylin - I - 1 - iyy.(2l + 1 + lyy. 



(48) 



with 



1/2 



Cni = -< — — — —3 } [(n + iy.\ ■ 

[ 2n[{n + iy.Y J 

(49) 

In the above-stated equations. En denotes the energy of 
the bound state to be studied, i.e., n = 2f; or n = 2e 
for the ground or excited states of the second electron, 
respectively. Since we are performing a qualitative anal- 
ysis, we will concentrate mostly on the functional form 
of Rni{r). The integral Ii present in the prefactor Vp^e 
can then be written as 

n^i (_i).+i2--i-i(^2^)-2-i r[2 + iy + 2l] 



v\{n - I - I - vy.{2l + I + vy. r[3/2-|-?] 



X2^^.(l + / + ^,^ + ^^+^,-M!). (50) 



The integral I2 in Ipie.kg is slightly more involved. It 
may be explicitly written as 

00 / 

^2-47r^ ^ {~i)'Yr{e,^,V,^)X2RX2n, (51) 



i=0 m=-l 



where 

1-2 R 

and 



rlRnjAr2)RnjAr2)3i{K.r2)dr2 (52) 



^20- j Yr{02,^2)[YC'{02,V2)YY;-^A02,^2)d^ 

(53) 

give the radial and angular dependencies of such prefac- 
tors, respectively. The explicit expression for T2 is then 



^2 = 4^ ^ ^ (-z)'(-i)™=rr(e«.,^K)i 



'{2lg + l){2le + l) 



47r(2/-f 1) 

'g, 'e, " I', "/ X'g, i-ej '"-g, — '"e "/ -^2H- (54) 

The radial integral T2R is proportional to 



/— /e|m— — / V 
X {lg,le, 0, |Z, 0) {Ig, le, mg, -nie \l, 0) I2R. 
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r [2 + A] 



lyeliygKne ^ le - I - Ve)\{ng ^ Ig - 1 - J/g)!(2;e + 1 + /^e)!(2Zg + 1 + J^g)! 



r[3/2 + /] VC" (^23,^^2 



1/2 



3+A 4+A 3 



2 '2 ' e{E2g,E2e) 



(55) 



where X — v^ + Vg + l^ + lg + l and C(-E'2g, ^'2e) is defined 
according to Eq. (f^. Note that the terms in Eq. ([5i)) 
are only non-vanishing if to = rUg — irie and li + I2 — I is 
even. 

In the present work, apart from the case in which only s 
states are involved and the angular integrals are constant, 
one may identify the following cases. First, the second 
electron may be initially in a p state and be excited to an 
s state. In this case, I = Ig = 1 and le = 0. Second, if the 



electron is initially in an s state and is excited to ap state, 
then Z = /e = 1 and Ig = 0. Finally, if the second electron 
suffers a transition from a p state to another p state, in 
principle / = 0,1,2. Due to the constraints upon / for 
the Clebsch-Gordan coefficients, however, only the terms 
with Z = 0, 2 will survive. Apart from that, the constraint 
upon TO will impose further restrictions for me and nig. 
The above-stated expressions, however, are applicable to 
generic hydrogenic states. 
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